Photon Splitting in a Very Strong Magnetic Field 
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Photon splitting in a very strong magnetic field is analyzed for energy u < 2m. The amplitude 
obtained on the base of operator-diagram technique is used. It is shown that in a magnetic field 
much higher than critical one the splitting amplitude is independent on the field. Our calculation is 
in a good agreement with previous results of Adler and in a strong contradiction with recent paper 
of Mentzel et al. 
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Virtual creation and annihilation of electron-positron 
pairs is known to induce nonlinear self-action of an elec- 
tromagnetic field. Photon splitting in an external field is 
one of corresponding processes of nonlinear QED. Obser- 
vation of a photon splitting is still a challenge for exper- 
iment. 

Theoretical study of this process has a rather long his- 
tory. Photon splitting in a constant and uniform exter- 
nal field was considered in the beginning of 70th in [1-4] , 
where earlier paper, containing errors, were cited. In 
101, the process was considered as a possible mech- 
anism for production of linearly polarized photons in a 
pulsar field (assuming the field H ~ Hq). At low pho- 
ton energies (a; <C m, m is the electron mass, we set 
h = c = 1) the splitting process can be analyzed by us- 
ing the Heisenberg-Euler (HE) effective Lagrangian. In 
the weak field limit {H ^ Hq), where Hq = m? /e = 
4.41 • IQ^^G is the critical magnetic field, the first term 
of expansion of HE effective Lagrangian can be used and 
the hexagon diagram contributes only. This was done in 
0, 1^. The polarization selection rules, especially with 
allowance for dispersion, were also obtained in see 
also textbook Sect. 129, 130 where the problem is 
given in detail. The comprehensive investigation of the 
process under consideration was carried out by Adler . 
For <C m and an arbitrary field strength the matrix 
element of the process was found as a result of the appli- 
cation of the full HE effective Lagrangian. At the same 
time, the allowed transition amplitude was calculated 
for the general case of an arbitrary field strength and 
photon energy below pair creation threshold {u> < 2m). 
A Green's function of the electron in an external mag- 
netic field in the Schwinger proper-time representation 
was used. Although the expression for this amplitude 
turned out to be very unwieldy for application, the am- 
plitude was calculated also numerically in wide interval 
of magnetic fields < H < Hq for = m as well as 
for uj <^ m. In photon splitting was considered in 
a crossed field E ± H., E — H, using likewise the elec- 
tron Green's function in the proper-time representation. 
Another form of the photon splitting amplitude in a mag- 
netic field in a general case was obtained in using sim- 
ilar approach for the Green's function calculation. Later, 
photon splitting in a constant and uniform electromag- 
netic field for arbitrary values of both field invariants 



was considered in [Q. The operator diagram technique 
developed by Katkov, Strakhovenko and one of us |^ was 
used. As a result, the solution of this technically quite 
cumbersome problem was substantially simplified. The 
amplitudes obtained in particular case of zero electric (or 
magnetic) field were found to be noticeably compact than 
those obtained in [||. The results of for a; <C m and 
H <^ Hq agree with those of [Q- jsj. In Q we performed 
numerical calculations for the case uj ^ m, because we 
were interested in another possibility: photon splitting 
in electric fields of single crystals at high energies ||] . It 
is worth to note that in all mentioned paper relativistic 
covariant and gauge invariant formulation of QED was 
used. 

Recently, photon splitting was considered once more 
]lO| . That was motivated by new astrophysics achieve- 
ments. In this calculation a non-covariant perturbation 
theory and Landau gauge were used. The results of this 
paper as well as its subsequent application are in a 
strong contradiction with all previous results. Matrix el- 
ement of photon splitting is found in JlOt in the form of 
very cumbersome threefold infinite sum. There is a series 
of short-comings in [ p^ : a) a low energy limit and weak 
field limit are not found from general expression; b) the 
authors of ]To| suppose that their approach is applicable 
for H = Hq, but the result of calculation with the full 
HE effective Lagrangian for a; <^ m is not reproduced; 
c) photon dispersion in a magnetic field is not taken into 
account; d) the study of cut off influence on summation 
with respect to Landau level numbers is not sufficient 
because the expression considered contains strong can- 
cellations. Adler |l^ criticized strongly papers [|lO| and 
pd| and suggested to make the independent calculation 
of photon splitting amplitude. 

Because of the potential astrophysics implications of 
the process (see e.g. |l^) we perform the numerical calcu- 
lation and analysis of the expression for photon splitting 
amplitude obtained using another formulation of QED in 
external field ||^ . We consider the most interesting region 
Lu < 2m and arbitrary H. 

Let photon with energy lu splits into two photons with 
energies lui and uj2. There is only one allowed transition 
in a magnetic field |]l| with respect to photon polariza- 
tions: B — > CC in notations of |7j or -L^|||| in notations 
of §. Putting electric field E = in general represen- 
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tation for photon splitting amplitude (eqs.(2.16)-(2.18), 0) we have for allowed transition: 
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t2{x — coshx — cosh(2i2 — 2:) 



, $ = [-2122(^1 ~ t2){tl -t2+x) - Ziti{ti - x)- 
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•22^2(^2 ~ a;)] /x — [(1 — 212:2) cosh a; — zi cosh(2ti — a;) — 22 cosh(2f2 — x) + 2122 cosh(a; + 2ti — 2t2)\ / (2 sinhx), 

G = [1 — (21 cosh(2ti — a;) + 22 cosh(2t2 ~ a;)) coshx] /x + 2cziZ2 sinh^(t2 — ii) [21 sinh^ ti + 22 sinh^(t2 — a;)] . 
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Here a — e'^ — 1/137, 2i_2 — ^1,2/^^, cr — sini?, 
is the angle between direction of magnetic field and mo- 
mentum of the initial photon. To derive (|l]) we rotated 
the contour of integration over each variable: x —>■ —ix, 
ti,2 —>- —iti,2- This transformation is valid for any 
Lj < 2m. As a result, the integrand in doesn't contain 
oscillating trigonometric functions. So, it is convenient 
for numerical calculation. It is easy to show that the am- 
plitude T is symmetric with respect to interchange of fi- 
nal photons (wi 0^2) ■ By virtue of gauge invariance the 
amplitude T oc wwiW2 if w — > and also T — > if — > 
for any u. That means that there are very strong com- 
pensations in and this circumstance should be taken 
into account under numerical integration. To overcome 
the difficulties it is convenient to perform subtraction in 
the integrand of (0): e^* e'^' 
G (proportional to 1/x) and e""" — > e' 
of the subtracted terms is equal to zero. 

For Lo m the main contribution to the amplitude is 
given by the domain of variables where c$ <C 1, C$0 1- 
Then expanding the corresponding exponents and keep- 
ing linear in c terms one can take the integrals over ti and 
^2. The result coincides with the photon splitting ampli- 
tude found with the use of full HE effective Lagrangian 
(eq.(22) in §). 

For u! ^ m the integrals in (P) are calculated numeri- 
cally. Without loss of generality one can put a = 1. In 
Fig.l the dependence of the amplitude T on the final pho- 
ton energy is shown for H = Hq/2 (a) and H — Hq (b) 
at different energies of the initial photon. The amplitude 
T is normalized on Tq: 
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For Lo/m = 0.1 the result coincides with a very good ac- 
curacy (better than lO^'^) with the result obtained from 
the full HE effective Lagrangian. 

The total probability of allowed transition vs w = uj/m 
is shown in Fig. 2, where 



Curve (1) corresponds to H/Hq — 1 and curve (2) to 
H/Hq = 1/2. Although the probability varies by many 
orders of magnitude in the interval of parameters con- 
sidered the essential part of the variation is absorbed by 
Wq. The function is nothing but the photon split- 
ting probability given by hexagon diagrams at w ^ m. 
Therefore, Fig. 2 shows the influence of higher order cor- 
rections with respect to uj/m and H/Hq. The probability 
found with the use of full HE effective Lagrangian also 
proportional to (uj/m)^. So, the intersection points of 
the curves with ordinate axis coincide with the proba- 
bility Whe found in this approximation. One can see 
from Fig. 2 that the probabilities W and Whe are essen- 
tially smaller than Wq at ~ Hq. At the same time, 
W/Whe grows appreciably at w — > 2m. Therefore, the 
exact photon energy-dependence should be taken into ac- 
count. Our numerical results agree (within a few percent) 
with obtained that by Adler in paper |^ where the case 
Lu = m was considered. 

The behavior of the amplitude T in a very strong mag- 
netic field H ^ Hq is of evident interest from theoretical 
point of view. In connection with this problem it is neces- 
sary to consider the selection rules for photon splitting in 
a strong field. Using the expression for eigenvalues of the 
photon polarization operator in a magnetic field found in 
Q, eqs.(3.33)-(3.35) we obtain that n|| = 1 -I- a/dir and 
nj_ =oc H for H 3> Hq and n± > n|| for any H. This 
means ( |^J^) that there is only one allowed transition, 
which we considered above, for any field H. 

For H Hq we found that the amplitude T is inde- 
pendent of magnetic field and can be evaluated in ana- 
lytic form. The main contribution to two-fold integral in 
(|l|) is given by the domain x ^ H/Hq and x — ^2 '^-^ 1- 
For threefold integral in (|l]) the main contribution comes 
from two domains: x ~ H/Hq, <i H/Hq, a: — t2 ~ 1; 
and X ~ H/Hq, ti ~ I, t2 ^ H/Hq. Performing the 
corresponding expansion, we obtain 
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where 

_ (47ra)-V^ 
^ 127r2 m2 

The aniphtude calculated using the full HE effective La- 
grangian a,i H is 

Tf/£; = TiWiW2/w^. The dependence of the amplitude T 
in this limit on the final photon energy (zi = loi/uj) is 
shown in Fig. 3 for different energies of the initial pho- 
ton. For strong field and uj — > 2m one can see in Figs.l 
and 3 a tendency of plato formation in the middle of the 
distribution. 

Thus, we performed the calculation of photon splitting 
amplitude using the exact formula valid for any magnetic 
field H and u) < 2m. If ^ to then our results coin- 
cide with the amplitude obtained from the full HE effec- 
tive Lagrangian. We obtained that in a very strong field 
H ^ Hq the amplitude doesn't depend on a magnetic 
field. We found that the refractive index n± > 7i|| for 
any H. So, there is only one allowed transition _L— >||||. 
Therefore, a photon cascade could develop only if mag- 
netic field changes it's direction (on distances much larger 
than the formation length of photon splitting). The re- 
sults of our calculation are in a good agreement with that 
obtained by Adler and in a strong contradiction with 
recent paper of Mentzel et al . 
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FIG. 1. The dependence of photon spUtting ampUtude on the final photon energy {zi = uji/uj) for H = Ho/2 (a) and 
H = Hq (b) at different energies of the initial photon: oj/m = 0.1(1), w/m = 1-5(2) and w/m = 1.9(3). The amplitude T is 
normalized on To given in the text. 
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FIG. 2. The dependence of total probability W (in term of Wo) on photon energy {w = uj/m) for H = Ho (curve 1) and 
H = Ho/2 (curve 2). The probability Wo is given in the text. 




FIG. 3. The dependence of the amplitude T on the final photon energy for H ^ Ho for different energies of the initial 
photon: oj/m = 0.1(1) , oj/m = 1.5(2) and oj/m = 1.99(3). 
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